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Abstract 

We study a geometrically frustrated triangular Ising antiferromagnet in an 
external magnetic field which is selectively diluted with nonmagnetic impu- 
rities employing an effective-field theory with correlations and Monte Carlo 
simulations. We focus on the frustration-relieving effects of such a selective 
dilution on the phase diagram and find that it can lead to rather intricate 
phase diagrams in the dilution-field parameters space. In particular, in a 
highly (weakly) diluted system the frustration is greatly (little) relieved and 
such a system is found to display only the second(first)-order phase tran- 
sitions at any field. On the other hand, for a wide interval of intermediate 
dilution values the transition remains second order at low fields but it changes 
to first order at higher fields and the system displays a tricritical behavior. 
The existence of the first-order transition in the region of intermediate dilu- 
tion and high fields is verified by Monte Carlo simulations. 

Keywords: Triangular Ising antiferromagnet, Frustration, Effective-field 
theory, Selective dilution, Phase transition, Tricritical point 

1. Introduction 

A triangular Ising antiferromagnet is a typical geometrically frustrated 
spin model which due to a high degree of frustration shows no long-range or- 
dering [l], [i[ . Nevertheless, a small perturbation, for example in the form of 
an applied external magnetic field or the injection of quenched non-magnetic 



"Corresponding author 

Email address: milan.zukovic@upjs.sk (M. Zukovic) 



Preprint submitted to Physica A 



December 24, 2012 



impurities, can lift the high degeneracy of the ground state and thus initiate 
long-range ordering. In particular, it has been shown that for a certain range 
of the field values the systems displays a transition between the ferrimagnetic 
phase with two sublattices aligned parallel and one antiparallel to the field at 
lower temperatures and the paramagnetic phase in which all spins are aligned 
parallel to the field at higher temperatures On the other hand, the in- 

jection of quenched magnetic vacancies in zero field can relieve frustration 
but only locally, which supposedly leads to a spin-glass order . Combin- 
ing both the presence of the external field and the magnetic impurities can 
lead to rather interesting phenomena, such as step-wise field- dependence of 
the magnetization or possibility of multiple reentrants in the phase diagram, 
as suggested by recent Monte Carlo [s| and effective-field theory [lOf studies. 
Another way of relieving the geometrical frustration is a selective dilution, 
in which different (in our case three) sublattices are diluted by non-magnetic 
impurities with different probabilities. Then the geometrical frustration can 
also be relieved globally, giving rise to long-range magnetic ordering phenom- 
ena even in the absence of the field. Kaya and Berker ll| have shown that if 
only one out of three sublattices of the triangular lattice is randomly diluted 
then a long-range order can develop in the remaining two sublattices already 
at relatively low concentrations of the vacancies. Various studies have shown 
that the presence of either the external magnetic field or the magnetic impu- 
rities has generally different effects on the critical behavior of non-frustrated 
and frustrated spin systems. For example, in the non-frustrated s pin sys- 
tems the critical temperature typically decreases with the field (e.g., [12H14J]) 



and in some peculiar systems, such as metamagnets [121 ] . the transition in a 
high-field and low-temperature region can change to a first-order one. On the 
other hand, in the frustrated systems the magnetic field can either enhance or 



suppress the transition temperature, depending on the field range [3H5|, [10 



Similarly, the introduction of the quenched dilution in the non-frustrated 
systems typically decreases the transition temperature down to zero at the 



percolation threshold [15[, while the opposite effect also known as "order by 
disorder" can be achieved by a selective quenched dilution in the frustrated 
systems [llj]. Therefore, by applying the external field to the selectively di- 
luted system we can anticipate different regimes of its effects on the critical 
properties of the system and thus an interesting phase diagram. 

The goal of the present study is to investigate changes in the critical be- 
havior of the system, the frustration of which is systematically controlled by 
the degree of the selective dilution between the fully frustrated (pure system 
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on triangular lattice) and non-frustrated (pure system on honeycomb lattice) 
limits. 



2. The model and method 

The selectively diluted Ising antiferromagnet in a field can be described 
by the Hamiltonian 

H = -jJ2 UASj ~hJ2 W 

where Si = ±1, are the Ising spins, h is the external magnetic field, J < is 
the exchange interaction parameter, and denotes the summation over 
all nearest neighbor (NN) pairs. & are quenched, uncorrelated random vari- 
ables which are equal to 1 with probability p when the site % is occupied by 
a magnetic atom and with probability 1 — p otherwise. Then p represents 
the mean concentration of magnetic atoms. For the current selectively di- 
luted case with only one (let us say A) sublattice diluted we will consider the 
sublattice-dependent concentrations < pa < 1, Pb = Pc — 1- 

In order to study the critical behavior of the system, we apply an effective 



field theory with correlations (EFT) (for review see, e.g., [16J). In contrast 
to the usual mean-field approach, EFT correctly accounts for the single-site 
kinematic relations through the Van der Waerden identity. As a result, EFT 
yields, for example, a non-zero critical concentration for quenched diluted 
systems, the lack of order in the one- dimensional Ising ferromagnet and the 
occurrence of order in the two-dimensional case with the critical tempera- 
ture improved over the usual mean-field theory |l6|]. However, an extra care 
must be taken when dealing with a frustrated system since a straightforward 
application of EFT may lead a complete loss of frustration and consequently 
meaningless results. Recently, we have applied EFT to an uniformly diluted 



triangular Ising antiferromagnet [10[ and in the present study we will follow 
the same procedure in order to include the geometrical frustration effects 
within EFT. Namely, we decompose the triangular lattice into three inter- 
penetrating sublattices A, B and C (see Fig. [TJ, such a way that spins on 
one sublattice can only interact with their NNs on the other two sublat- 
tices. Thus all the NN interactions are accounted for and the frustration 
results from the effort to simultaneously satisfy all the mutual antiferromag- 
netic intersublattice couplings. Such an approach has correctly reproduced 
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no long-range order behavior down to zero temperature for the pure system 
in zero field 10 . 



[Figure 1 about here.] 



2.1. Sublattice magnetizations 



The approach addopted in Ref. [lOj can be generalized by considering 



the concentration p sublattice-dependent and hence characterized by a three- 
component vector p = (pa,Pb,Pc)- Then, the A-sublattice magnetization per 
site can be calculated by taking a configurational average of the expression 



SW) = C{ n*i B i[£i cosh(/JJD) + £fSf sinh(/3 JD) + 1 - £?] 

x Ul A Mk cosh(/3JD) + sinh(/3 JD) + 1 - £ fc c ]) tanh(x + (3h)\ x=0 , 

(2) 

where Zab,2ac are the numbers of NNs of the spin from the sublattice 
A that belong to the sublattices B and C, respectively, (3 = l/k B T and 
D = d/dx is the differential operator. We can obtain similar expressions for 
the sublattices B and C from Eq. (J2]) by cyclic permutation of the indices 
A, B and C and for each sublattice taking an appropriate number of NNs 
on the remaining sublattices (in the present case zxy = 3, for all pairs X,Y 
G {A, B, C}). Then, we apply the decoupling approximation for the thermal 



averaging [10J and perform the configurational averaging of these equations 
for the respective sublattices. For the selectively diluted case considered in 
the present paper, i.e., (pa,Pb,Pc) = (pa, 1, 1), this leads to a system of 
coupled equations for the averaged sublattice magnetizations per site in the 
form 



|s=0) 



m A = Pa (a + bm-Qj (a + bmc) tanh [x + (3h) 
f^B = (oa + bm\) (a + bmc) tanh [x + (3h) \ x =o, (3) 
m c — i a A + brn^) 3 (a + bm^f tanh (x + (3h) 



=o- 



where a = cosh (f3JD), b = sinh ((3JD) and a a = 1 — Pa + Pa cosh ((3JD). 
By expanding the right-hand sides of Eqs. (jH]) we obtain 



4 



m A = K + Ki (m B + m c ) + K 2 (m| + 3m B m c + m 2 c ) + 

K 3 (m| + 9m B mc + 9memc + m^) + 
(mime + 3m|m^ + m B mc) + 

K 5 (m|mc + m|mc) + K^m B m%, 
m B = L + L x m A + L 2 m c + L 3 m A + L 4 m 2 c + L 5 m A m c + 

L 6 m A + L 7 m^ + L 8 m 2 A m c + L 9 m A m 2 c + . , 

Li m A m 2 c + Liim A m c + Li 2 m A mQ 

L Vi m A m 2 c + L 14 m 2 A m% + L 15 m A m^, 
m c = L + L\m A + L 2 m B + L 3 m A + L 4 m| + L 5 m A m B + 

L 6 m A + L 7 m B + L 8 m A m B + L 9 m A m|+ 

L 10 m A m B + L n m A m B + L 12 m A m B 

Li 3 m A m B + L u m 2 A m B + L 15 m A m B , 

where the coefficients i£j (z = — 6) and (j = — 15) are given in the 
appendix. 

Critical and tricritical points 

In order to determine the critical boundaries, let us introduce new vari- 
ables: the mean magnetization of the non-diluted sublattices B and C, 
nT>BC — { m B + m c) /2, and the order parameter mo = (m B — m c ) /2. The 
latter serves to distinguish between the ordered ferrimagnetic (mo > 0) and 
the paramagnetic (mo = 0) phases. Then, Eqs. (j3j) can be rewritten in terms 
of m A , m B c, and mo as follows: 



m A = A (m BC ) + A 2 (m BC ) m 2 Q + A 4 (m BC ) m G + A 6 m , 
m BC = B (m A ,m BC ) + 5 2 (m A ,m BC )m , (5) 
m Q = Ci(m A ,m BC )mo + C 3 (m A ,m BC )m , 



where 



(m B c) = + 2irim BC + 5K 2 m 2 BC + 20 J FT 3 m BC + 
5fQm BC + 2K 5 m BC + i^ 6 m BC , 
^2 (m BC ) = ~{K 2 + 12fT 3 m BC + 6fQm BC + 4FQm BC + 3K 6 m BC ) , (6) 
A 4 (m BC ) = K 4 + 2K 5 m BC + 3K 6 m BC , 
A% (m BC ) = -K 6 , 
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B (m A , wbc) = L + LirriA + L 2 m BC + L 3 m\ + L 4 m BC + L 5 m A m BC + 
L 6 m| + L 7 m BC + L 8 m 2 A m BC + L 9 m A m BC + 
L 10 m\m BC + L u m 3 A m BC + L 12 m A m BC 
L 13 m\m BC + L u m 2 A m BC + L 15 m A m BC , 

B 2 (m A , m BC ) = L 4 + 3L 7 m BC + L 9 m A + L 10 m 2 A + 3L 12 m A m BC 
Li 3 m A + 3L u m A m BC + 3L 15 m A m BC , 

(7) 

d (m A , m BC ) = ~{L 2 + 2L 4 m BC + L 5 m A + 3L 7 m BC + L 8 m 2 A + 

2L 9 m A m BC + 2L w m A m BC + Lnm A + 3L 12 m A m BC 
2L 13 m 3 A m BC + 3L u m 2 A m BC + 3Li 5 m A m BC ), 
C 3 (m A ) = -{L 7 + L 12 m A + L u m 2 A + L 15 m 3 A ) , 

(8) 

Now, by solving Eqs. §5$) in the neighborhood of a second-order transition 
line, where mo is small, and retaining the terms up to the third order of mo, 
one can obtain 

m A = m A0 +m Al m o + ..., (9) 
m BC = m B c + m BCl m 2 + . . . , (10) 
m = am + /3m Q + . . . , (11) 
where m Ao and m B c are solutions of 

m Ao = A (m B c ) , (12) 

m BCo = B (m Ao , m BCo ) , (13) 
and m Al , m B c 1 5 a an d P are given by 

m Al = E (m BCo )m BCl + A 2 (m B c ) , (14) 

Ei (m Ao , m BCo ) M (m B c ) + B 2 (m Ao , m BCo ) - 
mBCl " l-^(m BCo )^(m Ao ,m BCo )-E 2 ( mAo ,m BCo )' (15) 
a = d (m Ao ,m BCo ) , (16) 
(3 = C 3 (mA ) + E 3 (m Ao ,m BCo )m Al + E 4 (m Ao , m BCo )m BCl , (17) 

where 

cM (m BCo ) „ , N dB (m Ao} m BCo ) 

E (m BCo ) = — 1 , £i(m Ao ,m BCo ; - 



<im BC ' °' dm A 

c / x dB (m Ao ,m BCo ) . , <9Ci(m Ao , m BCo ) 

^ 2 (m A o,m B c ) = « , ^3(w Ao ,m BCo ; - 



<9m BC >l ' " ' " <9m A 



E 4 (m Aol m BCo ) = - 



dm B c 

The second-order phase transition line is then determined by a = 1 and 
(3 < 0. In the vicinity of the second-order phase transition line, the order 
parameter is given by 

m 2 = (1 - a)/0. (18) 

The right-hand side of the Eq. ffl8|) must be positive. If this is not the case, 
the transition is of the first order, and hence the point at which 

a = l,fi = (19) 



is the tricritical point (see, e.g., 17| and references therein). 



3. Results and discussion 

3.1. Ground-state configurations 

Let us first examine the ground-state spin configurations of the selectively 
diluted system with < p\ < 1 for different field values. Based on simple 
energy arguments it is easy to verify that, depending on the field value, there 
are three ground-state configurations. Their characteristics are summarized 
in Table [TJ The cases of pa = 1 and pa = correspond to the pure frustrated 
antiferromagnet on a triangular lattice and the pure non-frustrated antifer- 
romagnet on a honeycomb lattice, respectively (see Fig. [T]). The former 
three-sublattice system is known to display the ferrimagnetic phase with two 
sublattices aligned parallel and one antiparallel to the field within the field 



range < h/\J\ < 6 [3_H5|, [l(|, while the latter two-sublattice system shows 



an antiferromagnetic alignment within the field range < h/\ J\ < 3 [13l. Il4 
Therefore, for the field values h/\J\ G (0, 3], both these limiting cases display 
phase transitions, and we expect the phase transitions to occur also for in- 
termediate values of pa- For 3 < h/\J\ < 6, the system with pa = 0, i.e. the 
antiferromagnet on a honeycomb lattice, does not undergo a phase transition 



at any temperature [13l . Il4j . However, no phase transition can be expected 
even if pa > 0, since then the system will behave like a ferrimagnet. Namely, 
starting from the ground state spin arrangement (m A , tti-q, rnc) = (—Pa, 1, 1), 
the increasing thermal excitations will make spins on the diluted sublattice 
gradually flip in the field direction with no sharp phase transition. Appar- 
ently, for h/\J\ > 6 all spins are aligned to the field direction already at zero 
temperature, i.e. the system is in the paramagnetic state, and therefore no 



7 



phase transition can occur as the temperature is increased. For this reason, 
in the following we will focus on the critical behavior of the system in the 
region h/\J\ G [0, 3]. 

[Table 1 about here.] 

3.2. Critical and tricritical behavior 

In order to investigate the selective dilution effects on the critical behav- 
ior, it is useful to start with the fully diluted case of pa = 0, corresponding to 
the non-frustrated honeycomb lattice. Then, the increasing pa can be looked 
upon as the decoration of the honeycomb lattice by randomly placing the 
magnetic ions into the centers of the hexagons (Fig. [T]) and thus gradually 
increasing the frustration in the lattice. The effects on the critical behav- 
ior are apparent from Fig. [2j in which the phase boundaries are plotted in 
the field-temperature plane for different values of pa- For pa = 0, the sys- 
tem shows the phase transition within the field range < h/\J\ < 3, in 



accordance with the previous studies [13|, ]l4(, and the transition (or Neel) 



temperature at zero field knTiy/l Jl — 2.1038 coincides with the earlier EFT 



result for the honeycomb lattice [181 ] . With increasing pa the transition tem- 
peratures tend to decrease, however, no qualitative changes occur as long as 
Pa is relatively small, e.g., pa = 0.1. Nevertheless, upon further increase 
of pa, e.g., pa = 0.2 and larger, at sufficiently high fields the order of the 
transition changes from the second to first one at a tricritical point (TCP). 
Furthermore, for pa > 0.8 the character of the dependence of the critical tem- 
perature on the field changes from decreasing to increasing, and for pa ~ 0.8 
even both tendencies can be observed in the same curve within different 
ranges of the field. 

[Figure 2 about here.] 

[Figure 3 about here.] 

[Figure 4 about here.] 

Phase diagrams in the concentration-temperature parameter space are 
presented in Fig. [31 for several values of the field h/\J\. For h/\J\ = 0, 
a long-range order emerges below p\ = 0.8753 and the critical temperature 
increases with the decreasing value of pa up to kBT N /\ J\ = 2.1038 for pa = 0, 
corresponding to the honeycomb lattice. In a finite field, the phase transition 
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appears to take place for all concentrations, however, above a certain value of 
Pa the order of the transition changes to the first order at a TCP. With the 
increasing field the second-order transition temperatures decrease, as well as 
the concentrations pa corresponding to the TCPs. 

Hence, from the point of view of the varying degree of frustration, the 
following statements can be formulated. Relatively little frustrated systems 
with small values of pa < 0.1821 show only the second-order transitions. 
For more frustrated systems, within 0.1821 < pa < 0.9638, the transitions 
remain second order at low fields but they change at TCPs to first order for 
higher fields. Finally, for highly frustrated cases, with pa > 0.9638, there 
are no TCPs and the transitions are of first order at any field. In order to 
better visualize the behavior of both the tricritical temperature and tricritical 
field with the varying pa, we plot these dependencies in a separate Fig. HJ 
The tricritical behavior within 0.1821 < pa < 0.9638 is characterized by the 
monotonically decreasing tricritical field from Ht/\J\ = 3 at pa = 0.1821 
down to h T /\J\ = at pa = 0.9638. On the other hand, the tricritical 
temperature shows an inverted-U-shape dependence with the highest values 
of /csTr/| J | = 0.8682 in the middle and approaching zero at both ends of 
the interval 0.1821 < p A < 0.9638. 

3.3. First- order phase transitions 

Above, we established a line of tricritical temperatures below which the 
transitions are of first order. Even though the EFT formalism is able to 
determine the first-order character of the transition, unfortunately, it lacks 
an expression for the free energy and thus tools for its localizing, i.e., deter- 
mining the first-order transition boundary. Nevertheless, typical first-order 
transition features can be observed in the behavior of various quantities. In 
Fig. we plot temperature dependencies of the sublattice magnetizations, 
mA,mB,mc, the total magnetization, m = (niA + + mc)/3, and the 
order parameter, mo, in different regions of the parameter space. Appar- 
ently, there are qualitative difference among different cases. Indeed, while 



the curves in Fig. 5(b) , corresponding to a relatively low field and consider- 
able dilution, vanish to zero continuously, rather sharp drops typical for the 
first-order transitions are observed in the curves corresponding to either high 



fields (Figs. 5(c) and 5(d) ) or a relatively low field but no dilution (Fig. 5(a) ). 
More detailed computations in a narrow region around the point where the 
curves appear to be discontinuous reveal that the vertical drop results from 
the fact that above a certain value of the field the corresponding quantities 
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cease to be single-valued functions of the temperature. By exploring other 
solutions, we find the formation of wiggles, like those shown in Fig. E]for the 
order parameter, with h/\J\ = 2 and 2.5. Within effective-field approaches 



such a behavior signals appearance of a first-order transition [12l . [19 

[Figure 5 about here.] 

[Figure 6 about here.] 

At this point we note that our EFT results suggest that the pure sys- 
tem (pa — 1) shows first-order transitions for all fields. However, the pure 
system in a field was argued to belong to the same universality class as the 
three-state Potts model, which displays a second-order phase transition in 
two dimensions and the transition appeared second order also within 
the Monte Carlo mean- field approach j^]. Therefore, it is legitimate to ask 
whether the first-order transitions are not just an artifact of the EFT ap- 
proximation. Indeed, there are examples of other spin systems, such as some 
mixed-spin systems 2(J 21], for which EFT predicted the tricritical behavior 



but it was not confirmed by more accurate methods. 

[Figure 7 about here.] 

[Figure 8 about here.] 

In order to confirm the existence of the first-order transition in the present 
system, we additionally performed Monte Carlo simulations [22] . We note, 
however, that generally this is not an easy task since the first-order transition 
features often show up only when very large system sizes are used, which for 
highly frustrated systems means enormous equilibration times. For exam- 
ple, in our tests for pa = 1 and h/\J\ = 2.95 even as many as ten million 
Monte Carlo steps per spin (MCS) was not enough to reach equilibrium for 
the linear lattice sizes L ~ 200. Therefore, to demonstrate the existence of 
the first-order transition, we choose a case in which frustration is partially 
relieved, namely the system with = 0.5, and for which the EFT predicts a 
first-order transition (see Figs. [2]andH]). For this case, the relaxation times 
are reasonable and a fist-order transition character is observed at already as 
small linear lattice sizes as L = 12. In our MC simulations the number of 
MCS ranged from 2 x 10 4 , in the case of obtaining temperature dependences 
of various quantities, up to 10 7 , in the case of evaluating distributions at a 
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fixed temperature. First 10% of MCS were discarded for thermalization and 
the remaining were taken for thermal averaging. In Fig. [7] we show a tem- 
perature dependence of the order parameter for p& = 0.5 and h/\J\ = 2.95, 
for the increasing temperature. The curve displays a typical discontinuous 
first-order transition behavior with the jump above k B T j\ J\ ~ 0.6. A more 
precise (pseudo) critical temperature k B T N (L = 12)/|J| = 0.6545 is found 
by tunning the temperature so that the order parameter displays a bimodal 
distribution with equally high peaks (Fig. 8(a) ). The time evolution of the or- 
der parameter, shown for the first 10 6 MCS in Fig. 8(b), demonstrates that 
at fcsT/v(L = 12)/|J| the system switches between the ordered (mo > 0) 
and disordered (m Q = 0)0 states. The true critical temperature can be usu- 



ally estimated from a finite-size scaling (FSS) analysis 23j. Unfortunatelly, 



for frustrated spins systems increasing the lattice size would mean an enor- 
mous increase of the tunneling times needed to switch between the ordered 
and disordered states and thus hampering a reliable determination of the 
(pseudo) critical points needed for FSS by the equal height condition. We 
note that, although in our case we were able to observe clearly the bimodal 
distribution of the order parameter already at L = 12, relatively small lat- 
tice sizes generally suppress the typical first-order transition features. For 
example, when we used much larger lattice sizes, such as L = 240, some faint 
signs of the bimodal distribution appeared already for h/\J\ = 2.5. How- 
ever, as mentioned above, for such lattice sizes enormous equilibration and 
tunneling times would be required to reliably confirm the transition order. 
Thus, localizing a TCP, in which the order of the transition changes from the 
second to the first one, by the Monte Carlo technique would require using 
sufficiently large lattice sizes along with some more sophisticated method, 
such as the multicanonical Monte Carlo method, to keep the tunneling times 
in reasonable limits, and then performing careful finite-size scaling analy- 
sis 23[. Consequently, establishing the entire phase diagram, including the 
first-order transition boundaries and the tricritical points, in a broad field- 
concentration-temperature parameter space by Monte Carlo simulations is a 
computationally very demanding task and it is out of scope of the present 
study. 



1 Due to the fmiteness of the lattice size in MC simulations the order parameter remains 
non-zero even in the disordered phase and the zero value is only achieved for L — > oo. 
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4. Summary and conclusions 

We studied the critical behavior of a geometrically frustrated selectively 
site-diluted triangular lattice Ising antiferromagnet in a field within the 
framework of an effective-field theory with correlations. The selective di- 
lution was carried out by random removal of magnetic ions from one of the 
three sublattices. In particular, we focused on the effects of the frustration- 
relieving selective dilution on the resulting phase diagram. In accordance 
with some previous studies, in zero field such a controlled dilution relieved 
geometrical frustration and resulted in a long-range ordering in the remain- 
ing two sublattices below some threshold concentration of magnetic sites. 
In finite fields the system was found to display the phase transition at any 
concentration, however, the order of the transition depended on its value. 
In particular, for relatively little frustrated systems with small values of the 
concentration of magnetic sites we found only the second-order transition 
at any field. For more frustrated systems, corresponding to larger values of 
the concentration, the system displayed a tricritical behavior in which the 
transition remained second-order at low fields but it changed to first order at 
higher fields. Finally, for highly frustrated cases, with the concentrations ap- 
proaching the pure system, the transitions were identified as first order at any 
field. The existence of the first-order transitions in the region of moderate 
frustration and high fields was confirmed by Monte Carlo simulations. 

Nevertheless, the EFT approach only produced the second-order tran- 
sition boundaries and may have overestimated the region in the parameter 
space in which the first-order transitions are present. Therefore, addition- 
ally, it would be desirable to employ some other more reliable techniques to 
establish the first-order transition lines and also to check the extent of the 
parameter space in which they take place. 
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Appendix A. 

List of coefficients Ki (i = — 6) and Lj ( j = — 15) in Eqs. (j4]): 

iT = Pa° 6 tanh (x + (3h) \ x= o, 
K\ = 3pAa 5 6tanh (x + (3h) \ x =o, 
K 2 = 3pAd 4 b 2 tanh (x + (3h) \ x= o, 

K 3 = p A a 3 b 3 tanh (x + (3h) \ X=Q , (A.l) 
Ki = 3pA« 2 fr 4 tanh (x + (3h) \ x =o, 
K5 = Spxcib 5 tanh (x + f3h) | x= o, 
K 6 = p A b 6 tanh (x + (3h) \ x=0 , 

and 

-^0 = ai a3 tanh (x + (3h) \ x =o, 
Li = 3a A a 3 6tanh (x + (3h) | I=0 , 
L 2 = 3a A a 2 6tanh (x + | x=0 , 
L 3 = 3a A a 3 6 2 tanh (x + f3h) | s=0 , 
L 4 = 3a A afe 2 tanh (x + (3h) \ x =q, 
L 5 = 9a 2 A a 2 b 2 tanh (x + (5h) \ x=0 , 
L 6 = a 3 6 3 tanh (x + (3h) \ x=0 , 
L 7 = a 3 A b 3 tanh (x + (3h) \ x=0 , 

L 8 = 9a A a 2 6 3 tanh (x + (3h) \ x=0 , ^ ' ' 

Lg = 9a A a6 3 tanh (x + (3h) \ x =q, 

L w = 9a A a& 4 tanh (x + (3h) \ x= o, 

Lu = 3a 2 b 4 tanh (x + (3h) \ x =q, 

L12 = 3a 2 A b 4 tanh (x + (3h) | x=0 , 

L 13 = 3ab 5 tanh (x + f3h) \ x= o, 

L14 = 3a Ab 3 tanh (x + (3h) \ x =q, 

L15 = 6 6 tanh (x + /3/i) | x=0 . 

The coefficients can be calculated by applying a mathematical relation 
exp(XD)f(x) = f(x + X). 
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Sublattice: A- A B-Q C-O 

Figure 1: Triangular lattice partition into three sublattices A,B and C, where the sublattice 
A is randomly diluted. The signs + and — represent the spin values of +1 and — 1, 
respectively, and the question mark signifies the frustration arising from placing a magnetic 
A-ion in the center of the antiferromagnetically ordered hexagon of B- and C-ions. 
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Figure 3: Phase boundaries as functions of the concentration pa for different values of the 
field h/\J\. The circles represent the tricritical points. 
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Figure 4: Tricritical temperature and tricritical field as functions of the concentration pa- 
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Figure 5: Temperature variations of the sublattice magnetizations m\, and mc, the 
total magnetization m, and the order parameter mo, for different values of h/\J\ and pa- 
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Figure 7: Temperature variation of the order parameter mo for the increasing temperature, 
for p A = 0.5, h/\J\ = 2.95 and MCS= 2 x 10 4 . 
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(b) Time evolution of mo in (pseudo)critical temperature 

Figure 8: MC simulation results for the order parameter mo in the (pseudo) critical tem- 
perature UbTn /\ J\ — 0.6545, for L = 12, pa = 0.5 and h/\J\ = 2.95. Figure (a) shows the 
distribution of the entire history of 9 x 10 6 MCS and (b) shows the time evolution of the 
first 10 6 MCS, after equilibration over the initial 10 6 MCS. 
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Table 1: Ground-state sublattice magnetizations (mA,mB,mc) with the corresponding 
energies per site (H)/\J\N, for different field h/\J\ intervals. 



h 
\J\ 


(0,3] 


(3,6] 


(6, oo) 


(m A ,m B ,m c ) 


(PA, 1,-1) 


(-PA,1,1) 


(PA, 1, 1) 


W 

U\N 


1 hp A 
L 3 


1 2p A h[2 - pA} 


1 + 2p A h[2+ 3 PA) 
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